INT-PUB 04-01 



Dynamic universality class of the QCD critical point 

D. T. Son^'EI and M. A. Stephanov^'^-I 

^Institute for Nuclear Theory, University of Washington, Seattle, Washington 98195-1550 

^Department of Physics, University of Illinois, Chicago, Illinois 60607-7059 
^RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973 

(Dated: January 2004) 

Abstract 
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I. INTRODUCTION 



The phase diagram of QCD is a focus of many theoretical investigations Heavy- ion 
coUision experiments can probe a certain part of this phase diagram. The temperature and 
the baryon chemical potential at freezeout are determined, using statistical model fits to 
particle yields and spectra, to span the domain T ~ 100 — 180 MeV and /^b = 50 — 600 MeV 
This is the region where the critical point of QCD, as strongly suggested by model 
calculations and lattice Monte Carlo results |^, is located. The critical point is an end 
point of a first order phase transition line separating, in the chiral limit, the chirally broken 
and chirally symmetric phases. The precise location of the critical point is still unknown. 

Experimental signatures, based on the singular behavior of thermodynamic functions 
near the critical point, were suggested in Refs. |B, 0|. The characteristic feature of all 
such signatures is the non-monotonic dependence on the value of experimentally controlled 
parameter, such the critical point is approached and passed. An attempt to 

estimate the effects of critical dynamics on experimental observables was made in Ref. 0, 
where it has been shown that the effects can easily exceed the experimental background by 
orders of magnitude. The major limiting factor is the finite size and time effects, which 
round up the critical singularity in all observables. 

Near criticality, the crucial quantity is the value of the largest correlation length ^. For 
example, the singular contributions to event-by-event fluctuation observables studied in 
Ref. p are proportional to The divergence of ^ is limited by two effects: i) the finite 
system size, 0(10 fm) for heavy-ion collisions, and ii) the finite evolution time. As pointed 
out in Refs. 0|, the second effect is more important. Indeed, the time during which the 
the correlation length reaches its equilibrium value diverges as r ~ which defines the 
dynamic scaling exponent z. The finite evolution time limits the correlation length to be 
^ < (time)^''^. Since z > \ and typical evolution times of the heavy ion collisions are of the 
same order as the spatial size, the time limit is more stringent then spatial size limitation 
^ < (size). The value of z depends on the dynamic universality class of the critical point. 

The purpose of this paper is to determine this universality class. While the static uni- 
versality class is beyond doubt that of the 3d Ising model, the question about the dynamic 
universality class has not been satisfactorily answered in the literature. Some arguments 
were given in Ref. that the universality class could be that of model H in Hohenberg and 
Halperin's classification j^, - a model of the liquid-gas phase transition. The question was 
revisited in Ref. 0], where it was argued that the universality class is that of model C. The 
problem is not trivial because of an additional mode compared to the ordinary liquid-gas 
phase transition — the QCD chiral condensate qq. 

We shall see that mixing between the chiral condensate qq and the baryon density n 
leaves out only one truly hydrodynamic mode — a linear combination of the two. The hy- 
drodynamic mode is conserved and couples to the energy-momentum density. Consequently, 
the relevant hydrodynamic modes near the critical point are the same as in the model H. 
Thus we predict z ^ 3. 

Our results are in qualitative agreement with the arguments of Fujii based on a recent 
study p. The study P employed model calculations, and did not, in particular, take into 
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account the coupling to energy-momentum density. Our approach is mo del- independent and 
is based on the straightforward application of the hydrodynamic theory, reviewed, e.g., in 
Ref. Ei. 



The paper is organized as follows. In Sec. |n] we identify the set of hydrodynamic modes 
near the critical point. In Sec. II I II we discuss a simplified problem, where we forcibly ne- 
glect the energy-momentum density. The purpose is to demonstrate that only one linear 
combination of the chiral order parameter [qq) and the baryon number density is a true 
hydrodynamic mode. To keep the discussion simple, we shall work at the mean field level, 
sufficient to identify the correct universality class. We refer to the literature for going beyond 
mean field. In Sec. II VI we introduce the energy-momentum density, restoring the full set of 
hydrodynamic variables near the QCD critical point. We show that the relevant modes in 
the long-length, long-time limit are the same as in the universality class of the liquid-gas 
phase transition. We conclude with Sec. |3 In Appendix ^ we show that the inclusion of 
the isospin density does not affect the critical behavior. 



II. HYDRODYNAMIC MODES NEAR THE QCD CRITICAL POINT 

Our analysis is based on the premise that at sufficiently large distance and time scales, 
the dynamics of a finite-temperature system is described by a hydrodynamic theory. It is a 
theory which contains only the degrees of freedom varying slowly with time. Although the 
validity of hydrodynamics is amply supported by experiments, a generic theoretical deriva- 
tion, or proof, of the validity of such a description is a challenging problem. Nevertheless, 
the validity of hydrodynamic description has been verified in certain theories: scalar field 



theories jlJ| and strongly coupled gauge theories with gravity duals [1^ . We shall therefore 
assume that the dynamics of the finite-temperature QCD plasma can be described by a set 
of hydrodynamic equations. 

Generally, the full set of hydrodynamic modes include (i) densities of conserved charges, 
which relax via diffusion; (ii) phases of symmetry-breaking condensates; (iii) Abelian gauge 
fields of unbroken U(l) gauge symmetries. Near second-order phase transitions, (iv) the full 
order parameter (not only its phases, but also the magnitude) relax slowly, and therefore 
should be included in hydrodynamics. Thus, near the critical point of QCD on the (T, /i) 
plane, the modes potentially important for hydrodynamics are given by the fiuctuations of: 

• The conserved energy and momentum densities: e = — (T°°), and vr* = T*^*; 

• The conserved baryon number density, n = q'-f^q — {q'~f^q); 

• The chiral condensate a = qq — (qq) . 

The finiteness of the pion mass prevents pions from being a low-frequency hydrodynamic 
mode. This should be contrasted with the chiral limit, where the massless pions do become 
hydrodynamic modes, leading to a different universality class (both static and dynamic) . 

In QCD, the complete set of conserved charge densities also includes isospin density 
(neglecting small isospin-breaking effects). However, the coupling of the isospin to other 
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modes turns out to be irrelevant in the renormalization group sense (see Appendix IX]) . 
Therefore we do not need to include the isospin density in our treatment. 



III. WITHOUT ENERGY-MOMENTUM TENSOR 
A. Statics 

We start by discussing the statics. The static correlation functions can be found from a 
Ginzburg-Landau functional 

F[a,n] = j dx\^{d,af + hd,ad,n + ^{d,nf + V{a,n)\, (3.1) 

where 

A C 

V{a^ n) = —a'^ + Ban + —n^ + terms of higher orders. (3-2) 

The equilibrium distribution of the system is governed by e~^^ . The mixing between a and 
n is not forbidden at nonzero baryon chemical potential and quark masses, and has been 
taken into account. 

As we approach the phase transition, all the parameters of the Ginzburg-Landau func- 
tional (a, 6, c. A, B, C) remain finite, but the quadratic form in V becomes degenerate, 
AC = B^. At the phase transition, the potential has zero curvature along one direction in 
the (a, n) plane: 

n . a B C , , 

fiat direction: — = — - = — -— . (3.3) 

n A B ^ ^ 

The response of the system to static external perturbation can be found by including 
source terms into the Ginzburg-Landau free energy. For example, changing the chemical 
potential from to /i + and/or the quark mass from rriq to rriq + dm^ induces the 
following change in the average variables: 

a= A-\-B6i2- C (5mq ) , (3.4a) 

n= A-^{+ASfx+ BSm^). (3.4b) 

where 

A = AC-B\ (3.5) 

At the critical point A = and the responses are singular. In particular, the baryon 
susceptibility Xb = dn/dfi diverges at the phase transition. Near the phase transition, both 
changes induced by 6^ and Sniq occur along the flat direction ()3.3|1 of the potential. 

From the distribution e~^^ one can compute the fluctuations of a and n. In the limit 
q ^ the quadratic fluctuations are 

(^^-o) = ^ , (<^o) = TXB = '^. (3.6) 

Finally, by considering the response to nonuniform perturbations, q ^ 0, one can deter- 
mine the correlation length ^. Near the critical point it diverges as 
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B. Hydrodynamic equations 



At the linear order, the hydrodynamic equations can be written by using the standard 
rules reviewed in, e.g., Ref. l3] (see also Ref. We define the variables conjugate to cr 

and n: 

SF 

X„ = — = {A- aV^)(T + (5 - bV^)n , (3.8a) 

OCT 

5F 

Xn = ^ = (B - bV^)a + (C- cV^)n . (3.8b) 
on 

The linearized equations for a and n are then 

(^{q) = -laa{q)X„{q) - 7^„(q)X„(qr) + ^„{q) , (3.9a) 
n{q) = --fna{q)X„{q) - -fnn{q)Xn{q) + U{q) . (3.9b) 

We have written the equations for each spatial momentum q. Onsager's principle forces 
Ian = Ina- The uoise correlators are 

m,q)Ut\q')) = '^T^ijiq)i'^'^)'siq-q'W-t'), i,j = cT,n, (3.10) 

and are such so that the equilibrium distribution is given by e~^^ . 

In the limit of small momenta q ^ relevant to hydrodynamics, one can expand 'Jijiq) 
in powers of and keep only the leading terms. Due to the conservation of baryon charge, 
7„o- and 7„„ vanish in the limit q = and their expansions start at the order g^. There 
is no such constraint on 70-0-, so the expansion for this coefficient starts at the g° order. 
Introducing the notations 

7..(q) = r + 0(g2), (3.11a) 
laniq) = Xq' + Oiq^), (3.11b) 
Inniq) = + 0(g") , (3.11c) 

the hydrodynamic equations become 

a = -r^ + AV^^ + e., (3.12a) 
oa on 

n = AV^^ + AV^^ + a. (3.12b) 
da on 

with the noise correlators 

{U^)Uy)) = 2TT5\x - y) , (3.13a) 
{U^)Uy)) = -2T\6{t - t')V'6%x - y) , (3.13b) 
{Ux)Uy)) = -2TX6{t - t')V'6%x - y) . (3.13c) 
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C. Modes 



Now let us insert the expression ()3.H) for F into the hydrodynamic equations and find 
the dispersion relations for the normal modes. To leading order in the limit g — > we find 

a = -TAa-TBn, (3.14a) 
h = CXA + XB)V'^a + CXB + XC)V^n. (3.14b) 

The eigenfrequencies of the hydrodynamic equations are found by solving the equation 



det 



TA - iuj VB 
(XA + XB)q^ (XB + AC)g2 - lu 



(3.15) 



Near the critical point the two eigenfrequencies are 

uji = —iX—q^ . (3.16a) 
UJ2 = -iTA. (3.16b) 

Therefore, for small q there are two frequency scales: a small scale oc and a larger scale 
oc g". In the limit q ^ only Ui is truly hydrodynamic. It is a diffusive mode with the 
diffusion constant 

A 

D = X- = XC-X—. (3.17) 

A A ^ ^ 

The diffusion constant D tends to as the critical point is approached, A — > 0. To 
understand this fact consider A < 0. In this case the potential V has an unstable direc- 
tion. Negative D means that, instead of relaxing, perturbations grow, leading to spinodal 
decomposition. The term XB^/A in Eq. ()3.17|1 can be interpreted as "antidiffusion" — at 
the critical point it exactly cancels out the normal diffusion term AC 

To interpret the modes corresponding to lji and 002, we find the eigenmodes of the hy- 
drodynamic equations. Near q = the two eigenmodes are 

'a) (o)' 

At the critical point, the first mode is the fiat direction ()3.3|) of the potential energy V. 
Both cr and n vary along this direction. On the other hand, the mode with frequency UJ2 
corresponds to fiuctuations of a alone, unaccompanied by any change of n. 

The meaning of the two modes is simple. A generic smooth perturbation involving both 
a and n will relax over two distinct time scales. First, at a much shorter time scale, {rA)~^, 
the field cr{x) alone adjusts to the values which at each point x in space minimize the 
potential V{a,n) at a given local value of n{x): a = —{B/A)n. After that, at longer time 
scales, the chiral order parameter a can be forgotten. It simply "traces" the profile of n, 
which relaxes to ri = much slower, over the diffusive time scale {Dq'^)~^. From the point 
of view of slow dynamics at long distances, the addition of the chiral order parameter a 
does not change the physics, which is exactly that of the liquid-gas phase transition, with a 
conserved density as the only hydrodynamic mode (besides the energy-momentum). 
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D. Real-time correlators 



From the stochastic equations one can compute the real-time correlators of a and n. 
Assuming we are in the regime Dq^ <C TA, the correlators can be written as 

^ , ^ _ 2TT , 2TB'\q' 

Integrating over u we recover the equal-time correlators ()H.fi|) . Notice that while (cr^) obtains 
contribution from two characteristic scales u ~ TA and u ~ -Dq'^, the fluctuations of n are 
peaked only at the latter scale. 



E. Baryon diffusion, susceptibility and critical indices in d = 3 

We have seen that, from the point of view of real-time dynamics, the criticality is mani- 
fested by vanishing of the baryon diffusion rate D given by Eq. ()3.17p . Comparing to ()3.6p , 
we find that 

D = Axb'- (3.20) 

This relation has a simple physical explanation. Consider a configuration where n varies with 
space, n = n{x). This corresponds to a spatially varying chemical potential ijl{x) = XB^n^x). 
Recall that /i acts as the time component Aq of an external gauge field coupled to the 
baryon current. Thus such a chemical potential fj^{x) corresponds to an external electric 
field E = —'Vfi acting on the baryon charge. The baryon current induced by this field is 
jn = \E = — Axg^Vn, where A is the baryon conductivity. Using the definition of diffusion 
rate j„ = —Wn we obtain ()3.20|) . 

At d = 3 one has to take into account fluctuations which will modify the numerical values 
of the critical exponents. With the energy-momentum density fluctuations frozen, as in this 
section, the universality class is that of model B in the Hohenberg-Halperin classification 
It is a model of the uniaxial ferromagnet — with a single conserved mode. In the model 
B the diffusion constant is still inversely proportional to the susceptibility, as in Eq. ()3.20|) : 

^ ~ xb' ~ (3.21) 

but the divergence of xb is modified by r/ ^ 0.04 — a critical index of the 3d Ising model, 
relative to the mean field result. 

The dispersion relation u = —iDq^ applies only in the regime q ^ When q ^ 
we are in the critical regime, and higher powers of q are non- negligible. In this regime the 
scaling dictates dispersion uj ^ q^ . The two behaviors should match smoothly at g ~ 
From Eq. ()3.21|) one then finds z = A — r] (in mean-field theory z = A). 
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IV. COUPLING TO ENERGY-MOMENTUM 



So far we have completely neglected the motion of the plasma, regarding the latter as a 
static medium. Now we will allow the plasma to move — this will modify the hydrodynamic 
equations and change the value of the dynamic critical exponent. 

As we have learned in the previous Section, only one combination of a and n is truly 
hydrodynamic. This mode can be thought of as the baryon charge density n, with a simply 
tracing n and not independent. At sufficiently long timescales, the hydrodynamic theory 
contains this mode, the energy density e and the momentum density vr'. The linearized 
theory has four eigenmodes: the baryon diffusion mode, two diffusive transverse shear modes 
and a propagating longitudinal sound wave. 

At very low momenta, and long time scales, sound wave mode can be effectively integrated 
out, since it has linear dispersion ~ g and its frequency is much higher than the frequencies 
of the remaining diffusive modes ~ g^. In other words, the fluctuations of the energy 
density and pressure caused by sound excitations are fast and average out — there are no 
sound waves on these long diffusive time scales. The number of remaining hydrodynamic 
modes is three: two transverse components of vr* and the baryon density n} 

The same set of modes describes hydrodynamics of the liquid-gas phase transition jl5| . 
Here we onl y r eview the results of previous studies, referring the reader to the original 



literature [8|, |15[ for further details. 

The liquid-gas phase transition belongs to the dynamic universality class of model H. 
This model describes a system with a conserved order parameter, conserved and transverse 
momentum density, and nonzero Poisson bracket between the two. The same model describes 
the critical point of binary fluids. The hydrodynamic theory contains two kinetic coefficients: 
the shear viscosity f/ and the diffusion constant (in our case, that of the baryon charge D)? 

Both kinetic coefficients are singular (nonanalytic) at the phase transition: 

D = XXb ~ T^Xb', (4-la) 

n C""- (4.1b) 

The two new dynamical critical exponents satisfy the following relation 

XA + x^ = 4-rf-?7. (4.2) 

The enhancement of D by a power of is due to the contribution of convection to the 
baryon conductivity A. Although the calculation of x\ and individually is not simple, 
the relation ()4.2j) has a simple physical explanation Let us apply the field E = — V/x 



^ The physically intuitive argument presented here is corroborated by the analysis of linear mixing of four 
variables: cr, n, e and the longitudinal component of tt*. As in Section lTTTl there is only one diffusive mode, 
in which a and e trace the local value of n to minimize the local value of the thermodynamic potential 
(i.e., maximize the pressure). 

^ We choose the notation fj for the shear viscosity to avoid confusion with the static critical exponent rj. 
The constant of heat conductance k is not an independent constant, but is related to I? by Z? oc k,Xb^ ■ 



8 



introduced already in Section UTTl In model B, without convection, this field would induce 
baryon current via diffusion: j„ = XE, where A is finite. However, in this field E, the 
baryon charge carrying fluid will experience mechanical force equal to nE per unit volume. 
The fluid will accelerate to velocity v at which viscous drag balances the external force: 
/vise + /appi = . For a chunk of fluid of typical linear dimension L the drag and the apphed 
forces are of order: 

Use -VPvL''-^ and f^^^.^nEL'', (4.3) 

where p is the mass density. The corresponding velocity v diverges with L and so does the 
induced baryon current: 

jn = nv r-. —L^E. (4.4) 
VP 

The divergence is cut off at the scale L ^ because the baryon number fluctuation, given 
by (n^) = Tx-q/L'^ for L ^ ^, is correlated on at most that scale. Therefore at the critical 
point we flnd a singular contribution to the product of the kinetic coefficients 

\n ~ {n^)e - Xb^-" - (4.5) 

which leads to Eq. 

In the (i = 3 Ising model r/ is very small. It has been found BHQ that the divergence 
of the shear viscosity is very weak, or is numerically small. According to Eq. ()4.2p . 
therefore, a; ~ 1. Calculations using renormalization group in 4 — e dimensions yields the 
values 



)ie, XX 



= — e(l + 0.238e + ■■■)- 0.065 , (4.6a) 
19 

XX = — e(l - 0.033e + ■■■)- 0.916 , (4.6b) 
19 

while mode-coupling calculations done directly in = 3 and assuming 77 = yield [l^ 

= 0.054, xa = 0.946. (4.7) 

One sees that in model H the diffusion constant still vanishes at the critical point: 

D ~ ^-^+^+^\ (4.8) 

but the power is closer to 1 rather than to 2 as in model B. That is due, as we have mentioned, 
to fluctuations of the fluid motion. Using the same argument as the one presented at the 
end of Sec. 11111 we flnd 

z = 4:- rj - xx ^ S. (4.9) 



V. CONCLUSION 

We have seen that the dynamic universality class of the QCD critical point is that of model 
H, i.e., the liquid-gas phase transition. The chiral order parameter, being non-conserved and 
mixing with the conserved baryon charge, does not affect the dynamic universality class. 
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Here we would also like to compare our results to the universality class argument of 
Ref. (7|]. The argument is based on three assumptions, spelled in Ref. i) the chiral order 
parameter a is not conserved; ii) there are other conserved quantities, such as the baryon 
density n; iii) the Poisson brackets between a and the conserved quantities vanish. While 
the first two assumptions are correct, the last assumption is not. There is a non-zero Poisson 
bracket between a (or n) and the momentum density [7ii{x),o'{y)] = a{x)'Vi6^{x — y). As 
we have seen, the coupling to momentum density plays essential role in determining the 
dynamic universality class. Another, more subtle implicit assumption of Ref. 01 is that the 
order parameter a and baryon density n cannot mix. This assumption leads to model C, 
where the mixing is forbidden by a symmetry of the order parameter (it would be a — >■ —a 
in this case). As we have seen in Section ITTT| the an mixing eliminates nonconserved mode 
from hydrodynamic theory, leading instead to model B. 

The model H value of z ~ 3 is larger than the model C value z ^ 2.17 used in Ref. 0|. 
As discussed in the introduction, this means that the effect of the time constraint on the 
correlation length ^ is stronger. Thus the numerical estimate of maximal ^ in Ref. should 
be revised downward. However, we do not expect the main conclusions of Ref. to change 
qualitatively.'^ 

An interesting step beyond the relaxation equation for the diverging correlation length 
in Ref. 0, which one can contemplate, having the correct hydrodynamic theory in hand, 
is to consider real-time evolution of the hydrodynamic modes near criticality, similar to the 
study of Ref. 0. 

The pions do not enter the hydrodynamic theory as long as the correlation length ^ is 
sufficiently large compared to the inverse of the pion mass 1/m^ (see Section |n]). Although 
this condition is always fulfilled sufficiently close to the critical point jin the realistic case of a 
heavy ion collision, the maximal achievable ^ (2 — 3 fm, according to ja,l2|) is not significantly 
larger than l/m^r. Therefore it might be interesting to study the effects due to the crossover 
between the two regimes ^ ^ ^/ttIt^ (model H) and ^ <^ ^/miT (0(4) antiferromagnet (l3|). 

It would be interesting to explore the phenomenological consequences of the vanishing 
baryon diffusion rate for heavy ion collisions. Most likely it will manifest itself in the 
fluctuations of the baryon number jl9| . 
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APPENDIX A: ISOSPIN DENSITY AND SCALING DIMENSION COUNTING 

In this appendix we show that the inclusion of the isospin density rij does not affect the 
critical behavior. For simplicity, let us neglect the energy-momentum tensor and discuss only 
the coupled system of the isospin density and the baryon density. The coupled dynamics at 
the critical point is given by a pair of stochastic equations: 

fiB = Ab (-cV^^B + gBnl + ^^i) + ; (Ala) 
hi = Ai (xi^ni + ginf + gBiiriBni^ + 6 • (Alb) 

Note that the term AbXb^^^^b in Ecl (lAlaj) vanishes at the critical point, because xb — ^ cxd. 
The isospin susceptibility xi is finite jl^] . We kept only couphng terms of lowest order (most 
relevant in the infrared) consistent with the symmetry rij —>■ —ni. The noises ^b and are 
auto-correlated as 

iUt, x)Ut', y)) = -2T\B5{t - t')\/'S''ix - y) ■ (A2a) 
(ei(t, x)Ut\ y)) = -2T\i6{t - t')VH\x - y) . (A2b) 

Here d is the number of spatial dimensions. 

The canonical scaling dimensions involved in renormalization group transformations are 
determined as follows. Comparing linear terms in Eq. ()Ala|) we establish that if space has 
dimension —1, [x] = —1, then the dimension of time is [t] = —4. The dimension of noises 
are then found from Eqs. ()A2j) to be [^b] = [6] = d/2 + 3. From Eq. (jAlajl one then finds 
[hb] = d/2 — 1. In Eq. (jAlbj) we see that the term hi can be neglected in the infrared, so 
[m] = d/2 + l. 

Using these scaling dimensions for the fields, we can now determine the scaling dimensions 
of the coupling terms. One finds that [ABf^Bii] = —1 — d/2, and [Ajfy'Bii] = I — d/2, so 
the coupling between baryon and isospin densities is irrelevant for d > 2. (In contrast, 
[AbS'b] = 4 — d, so the n| term is relevant for d < 4.) 

Similarly one can show that isospin density remains decoupled from the critical dynamics 
when the energy-momentum tensor is included. This simple power counting scheme can be 
formalized by rewriting the stochastic equations as a path integral. 
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